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Abstract 

We consider the thawing model in the framework of coupled quintessence model. The effective 
potential has symmetry which is broken spontaneously when the dark matter density becomes 
less than a critical value leading the quintessence equation of state parameter to deviate from -1. 
Conditions required for this procedure are obtained and analytical solution for the equation of 
state parameter is derived. 


1 Introduction 

Current observational evidences indicate that our universe, consists of nearly 68% dark energy, the 
mysterious energy component with negative pressure, and nearly 27% of an unusual matter component 
called dark matter [T] . Although , considering dark energy as the cosmological constant A was proposed 
in standard A cold dark matter, ACVM model [2], due to problems such as the coincidence and 
fine tuning, it may be more interesting to think about dark energy as a dynamical component. For 
this purpose, dynamical scalar fields in the framework of quintessence models are suitable 

candidates, which have been widely studied before [3]. In ACVAi model, the equation of state (EoS) 
parameter of dark energy is exactly —1 while observations show the EoS parameter of dark energy, 
W(f,, satisfies —0.9 < < —1.1 [1], which means that may not be exactly -1. One possible way to 

explain this is to consider dark energy as “Thawing quintessence” [1]. Thawing scalar fields hardly 
move in the past but roll down the potential recently. These type of models are studied in detail in [5] 
but without considering any interaction between dark sectors. To alleviate the coincidence problem, 
various interactions between dark matter and dark energy have been proposed in the literature [6]. In 
our study we consider coupled quintessence thawing model. We use this coupling to explain deviation 
of wp, from —1 through spontaneous symmetry breaking of the effective potential. 

In our investigation we consider a coupled dark energy model in which the effective potential of 
the scalar field depends on dark matter density such that the scalar field resides in the minimum of 
the effective potential during the time when the density of matter is greater than a critical value [7]. 
However, when the density becomes less than a critical density, the shape of the potential changes 
and, as we shall see, the initial symmetry spontaneously breaks and the scalar field rolls down to a 
new minimum and the equation of state of (j) deviates from -1. Unlike E we do not restrict ourselves 
to a specific potential or coupling and do not consider any additional coupling between dark energy 
and baryonic matter. So we relax the local gravitational tests and the gravitational screening effect 
which pose a large mass constraint on the quintessence mass. As the nature of dark sectors is not 
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specified, considering such phenomenological interactions between (only) dark matter and dark energy 
is customary in the literature. 

Our main goal is to derive an analytic expression for the equation of state parameter of the 
quintessence in terms of the scale factor near the maximum of the potential and study analytically 
the deviation of w^j, from — 1 after the symmetry breaking. This is an extension of studies on thawing 
models in the literature where such an interaction was not considered [5] and the reason behind the 
scalar field leaving its initial position was not clearly shown [ 8 ]. 

The plan of the paper is as follows: In the second section, we discuss some preliminaries regard¬ 
ing our coupled quintessence. In the third section, to establish our model, we propose the required 
propositions which are necessary to have the desired symmetry breaking. In the Sec.4, which contains 
our main results, we study analytically the evolution of the scalar field and its EoS parameter near 
the local extremum of the potential where ~ — 1 and plot the EoS parameter of dark energy 
numerically. We use units h = c = 8 ttG = 1 through the paper. 

2 Preliminaries 

We consider a universe described by a spatially flat Friedmann-Lematre-Robertson-Walker (FLRW) 
metric, filled by a quintessence scalar field dark energy (f, cold (pressureless) dark matter p^, and 
baryonic matter p^. The dark energy is coupled to dark matter via a source term, f{(l))cj)pd [9]. For 
example One can see such a coupling after a conformal transformation of Brans-Dicke theory or in 
string theory m- We consider the interaction only between dark sectors and the baryonic matter 
does not couple to dark sectors 


^ + 3H<j) + V,j> = -f{(t))pd 
Pd + 3Hpd = fif))<j)pd 

Pb + 3Hpb = 0. (1) 

The Friedmann equations are 

H^ = l(^l<p^ + V{cf) + Pb + Pdy (2) 

and 

H + Pb + Pdj ■ (3) 

Note that dU), and only one of the two equations ([2|) and ([3]) can be considered as independent 
equations. H = ^\s the Hubble parameter and a{t) is the scale factor. By convention, we take a = 1 
at the present time. The equations of motion can also be obtained via the action m 

5 = / d‘^x^/^ (^^R-^g^^’'d^(pd^4>-V{(f)^ + Sm[g\u^'fp"], (4) 


where R is the scalar curvature, (p is the quintessence and Sm is the action corresponding to other 
species if'' . The quintessence is coupled to via g^^^ = A‘^{(f>)g^y, where Ai{cj)) is a positive function. 
We take different couplings for different species m- As we aim to consider interaction only between 
dark sectors, we assume that Aj is nontrivial Aj 7 ^ 1 only for dark species. If for dark matter we take 

Ai = A{(p) = exp f{4>)d(pj , (5) 

by variation of the action we get the equations of motion © and Friedmann equation ([2]). 

By introducing the density, p, 

Pd = A{(t))p, ( 6 ) 
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we find that 


P + 3-H'/9 = 0, (7) 

whose solution is given simply by /? = where p{a = 1) = p^. This p is not the physical 

density but a mathematical object which help us to solve and interpret some of our equations. The 
quintessence equation of motion can be rewritten as 

(j) + 3H^ + = —A^^p. ( 8 ) 

Hence in the evolution of (f>, V^^^{(f)^p) defined by + A^^p, plays the role of an effective 

action. 

ij) + 3H^ + = 0. (9) 

It is worth noting that V^^^{(l),p) depends on (j) as well as on p. Therefore the value of f which 
extremizes the effective potential, is not a constant generally and may be p (and consequently time) 
dependent. 

3 Conditions required for thawing via symmetry breaking 

In this part, necessary conditions for A((j)) and V{(j)) allowing the thawing process via symmetry 
breaking are investigated. Usually in the literature, some specific examples for V{(f)) and j 4((^) are 
used to study the symmetry breaking in symmetron models [7], but we try to discuss the subject 
generally. 

To attribute the evolution of w{4>) in thawing model to spontaneous symmetry breaking in our 
coupled quintessence model, following [7j, we assume that the effective potential has a Z 2 symmetry 
which is broken by the evolution of p. We consider two following eras: 

I- When p is greater than a critical value, the quintessence field stays at the minimum of its 
effective potential. 


II- When p becomes less than a critical value (specified by the parameters of the model)the shape 
of effective potential changes and the Z 2 symmetry is spontaneously broken. The previous minimum 
becomes an unstable point, and the scalar field moves towards the new minimum of the effective 
potential. 

To study the condition (I), let us investigate whether the scalar field can reside at the minimum 
of the effective potential {cf), p{f)), in the sense that 

U</) + -^,<pP ~ 0- (10) 

As p = poa~^, let us see whether (fTOjl gives a constant or time dependent solution for (p. Note that 
since we consider m as an exact equation, in contrast to works where (HOD is only considered as an 
approximation, using the equation of motion ([ 8 D, we immediately find 

p + 3Hf) = 0. (II) 

We consider two possible situations (i and ii) and show that the first one (i.e. i) is not acceptable, 
i- A^,p 7 ^ 0 at the minimum of the effective potential: In this situation (jlOp determines p in terms of 

4 >, 

P = (12) 
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where h{4>) is an analytical function of (j). The solution of dill) is 

^ = Ca-^ 


where C IS mimeriCcil constctiit. R,ec 3 /llirig thcit p ssitisfies (j7p ctiid by usiii^ (jl2p cnid m , we can 
obtain an integral equation for the scale factor 


h{ j Ca-^dt) 


— Poo. 
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(14) 


But this is not the whole story: We have one more equation to he satisfied, i.e. the Friedmann equation 
m or m- As the Friedmann equation is independent of the equations ([7]), ([2]), it is natural to expect 
that it may be not consistent with dm)- To clarify this issue, let us consider, as an example, a 
quintessence held with the potential 


+ v, 

where A and v are two positive constants , and 

Note that p = h^fpcj). Using ([7|), this leads to —2>Hp = h^^pcj), whose time derivative gives 

- 3Hp - 3Hp = + h^^^. 

But from p = h^^pcj) we have —3Hp = —which by using m gives 

— 3Hp = h^^cp. 

Now, using (USD, (HID reduces to 


— 3Hp = h 
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(15) 


(16) 


(17) 


(18) 


(19) 


For (USD and (USD, /?-,</,</> = [-m — ~ 


M2 


= < 0, therefore from (flSj) we hnd that H > 0. This 


is in contradiction with the equation of motion ([ 2 D 


H = — 
2 


+ Pb + 


( 20 ) 


which dictates that H is negative in our model (note that the quintessence is assumed to be real and 
A{(p) > 0). Despite this, note that in the presence of interaction it is possible to obtain an effective 
equation of state parameter for the dark sector such that < — 1 [IS]. In the same way, it is not 

hard to see that the assumption (i) is not consistent with the Chameleon model m- 


ii- = 0. However if we assume that at the minimum of the effective potential, = 0, then (|10l) 
does not put any constraint on p nor on a (note that p / 0 ). As A is a specific function of (p, = 0 

fixes the value of <p, i.e (p = (pc, where (pc is a time independent constant. Hence C = 0 in (fT^ . Besides, 
from dinD, we conclude V,f){(pc) = 0. In this case, the potential of the scalar field plays the role of 
a cosmological constant and the Friedmann equation holds consistently. For example, for the model 
specihed by (USD and dUD, the quintessence may lie at (p = (pc = 0, initially, where A^ 0 (O) = 0 , but 
for Chameleon model with A{<p) = exp{/3<p), the quintessence can not reside at the minimum of the 
effective potential in any situations. 

Whence at the stage (I)(mentioned in the beginning of the section), (p must be a constant: (p = (pc, 
<p = 0, and A^^{(pc) = V(f,{(pc) = 0. In this era the EoS parameter of the quintessence is 




h4>^ - V{f>) 
W + V{(P) 


( 21 ) 
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as required. 

It is important to note that (I), which states that the sealar field resides exaetly at the the minimum 
of the effeetive potential, is mueh more restrictive than the adiabatieity situation used in the literature 
for models such as Chameleon, where the field traces the minimum of the effeetive potential m 
and m is only an approximation and so m does not hold and higher derivatives of the effective 
potential should he considered. 

In the first stage, i.e. stage I, as = 0, then ^ = 0 has a constant solution 

(j) = 4>c- In this stage pd and pb evolve similarly; pd (x pb oa a~^. To see whether the system is stable in 
the (p direction, i.e. small quantum fluctuations do not force (p to leave pcj it is sufficient to consider 
second derivative of the effective potential. As long as + A,(j)(j)P) 1^^^ > 0 (the shape of 

the effective potential is convex), p continues to stay classically at p = pc- In this stage the vacuum 
expectation value of the scalar field is {p) = pc- But during this epoch p decreases, and if we establish 
our model such that V,p(j){pc) < 0 and A^^,p{pc) > 0, we eventually have < 0. At this 

moment, the scalar field becomes tachyonic and pc is no more the true vacuum, and p recedes from 
pc and moves towards the true vacuum, pd- Finally we will have < p >= pd after a time of order 
Tufj-j,. In this stage p 0 and the EoS parameter deviates from w,j, = —1, due to the mentioned 
symmetry breaking. 

To see, via a dynamical system analysis which includes the late time behavior of all ingredients 
of the system, that the quintessence cannot stay at the initial point forever (playing the role of a 
cosmological constant), one can also use a simple and brief phase space analysis. Note that in this 
brief discussion we only consider the late time evolution, and the (transient) stability of the scalar 
field which occurs at the first stage (I) and studied above, is not discussed in this dynamical analysis. 

By defining dimensionless variables [H] x = U = 


dx 

dN 

dzd 

dN 

dzb 

dN 

dN 

dg 

dN 


-3a; - + 3x3 3^2^ 3^2^ _ ^zjg, 

r- 2 3 3 o 3 2 

VGxzdg + SzdX - -Zd + -z^ + -ZdZ^, , 

2 3 3 o 3 2 

3x Zb - -Zb + -z^ + -ZbZ^, 

V6xf (T - 1), 

V6xg^ (C - 1), 


( 22 ) 


where = Ina, / = 


_ V' 


4 ‘ 


= T = 

- /t •! -L — 


vv" 


C = 


_ AA" 
-J!T 


From 


we have x^ + + z? + z? = 1, 


jr, y — ^ — -yJT, S = ~A^- nave a. -r y -r ^d 

SO y is not considered as an independent variable and may be rewritten in terms of other variables. 
Therefore (j22p becomes an autonomous system of differential equations by taking F = T{f,g) and 
C{f,g). The point D = {x = 0, Zd = 0, Zb = 0 }, / = f{p) = 0 is a critical point of the system, where 
the energy of the Universe is coming only from the scalar field potential, at late times. This point 
is exactly the same point of the phase space if the scalar field could stay at p = pc- Note that the 
relations Zd = 0, Zb = 0 hold only at late times and are not true in the first stage where > 0 . 

However this can happen provided that D is a stable point. To study this stability let us consider 

f = f + df. 

Consider matrix M defined by 


A 


X + 6x, 

Zd = Zd + dzd, Zb 

/ <5x \ 


/ dx \ 

dZd 


dzd 

6 Zb 

= M 

dzb 

df 


df 

\ dg y 


\ dg / 


(23) 
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Hence at the critical point D 


M = 


-3 
0 
0 

^ V{<P) 

V VGf{C{0,9)-l) 


0 

0 

0 


0 

0 


0 

0 


o\ 


0 

0 

0 


0 0 / 


(24) 


Eigenvalues of M are A = 0, A = —A = —A = —| ± \ • Assuming V{4>) > 0, which 

is necessary to have a real Hubble parameter, one obtains a positive eigenvalue when V"{<p) < 0, 
indicating the instability at the point D where the potential is concave (note that at point D where 
Pd and pb tend to zero, the effective potential becomes asymptotically the same as the potential at 
late time). 

Note also that condition V"{^) < 0, which leads to a positive eigenvalue i.e. A = — > 

0, is a sufficient condition for having instability at cj) = (pc at the second stage. As we require that at 
late time the scalar field becomes unstable at p = pc, we adopt V”{p) < 0 for our proposed model. 
However this is not a necessary condition for occurrence of instability; as one of the eigenvalues is 
zero, the critical point is non-hyperbolic and even if the other eigenvalues are negative, one must 
apply other tests such as applying center manifold theorem to check the stability or instability of the 
system at the critical point m- So it may be possible that the critical point is unstable at late times 
even for V"{p) > 0. 

In summary, to satisfy our proposal H, we assume that the scalar field stays at the maximum of 
its potential (where V = 0 , V < 0 , and > 0 ) initially. This is the minimum of the effective 

potential and this stay is transient. Later, when this point is no longer the minimum of the effective 
potential, due to instabilities, the scalar field moves and we obtain a dynamical held, that \s p ^ t). 
From m, it is clear that > — 1 . 

Theories like (H]), with coupling between dark matter and the quintessence may exhibit instability 
in the adiabatic regime where the quintessence follows the minimum of its effective potential |16j . The 
sound speed squared becomes negative and perturbations grow exponentially. In by considering 
the action Q, it was shown that such instability occurs only for coupling satisfying 


dln{A{p)) J_ 

dp Mp' 


(25) 


where Mp is the reduced planck mass (we have taken Mp = 1 in the units used in the paper). 
So one can evade the instability by choosing appropriate parameters. For example if one chooses 
A{p) = 1 + ^ + O , where M is a mass scale such that ^ <C 1, the condition (1251) reduces to 

p S> ^. So for p < ^ we can evade this instability. 

Summarising, the model ([T]) can describe the thawing model via Z 2 symmetry breaking. For this 
purpose, it is sufficient to take an even potential and even A{p) which at the initial point pc satisfy 
y,(l){Pc) = 0 , A^^{pc) = 0 , V^^ipc) < 0 , and A^^^(pc) > 0. 


4 Evolution of the model near the extremum of the potential 

In this section we begin by deriving an analytical expression for the EoS parameter of coupled dark 
energy. Although p plays the role of a cosmological constant initially, afterwards, by the symmetry 
breaking, the scalar field rolls down from pc and deviates from —1. The difference of this model and 
ACTyjti is that it describes a dynamical dark energy with w > —1 while for ACDJti we have always 
w = —1. This may help to alleviate the coincidence problem and also improve the compatibility of 
the model with astrophysical data m- 
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In the following we study evolutions of the scalar field (p and near (p = (pc, when ~ — 1. 
We assume that the present epoch is in this regime and take oq = 1. By the subscript 0, we denote 
the value of the corresponding parameter at present time. Note that our computation is only valid 
for the present epoch when assume that the scalar field is near pc, that is 0 ~ (/>c. In the absence of 
interaction, one obtains [18] 

= Hi ( 11^0 + (1 - . (26) 

Relative densities are defined by Oj = e. g. = ^ 7 , where p^j, is the energy density of the 
scalar field. As ~ — 1 , we have ~ 0 implying that ~ ^ 4 , 0 , which was used in derivation of 
(l26l) . The solution of (126]) is [18] 


a = 


/1 - \ 

V ^(/>0 / 



(27) 


where t\ = 




In our coupled quintessence model we have an interaction term between dark sectors 


Pd + ?>Hpd = A ^ 4 , A ^pdy/l4>P4„ 

P4) T ^H'^ 4 )Pfp A^ 4 )A Pd's/'^pPp, ( 2 S) 

where ' 44 , = W 4 ) + 1 and we have used p"^ = '^ 4 >P 4 )- Note that in the regime under study 7 ^ ~ 0. The 

_3 

scalar field satisfies equation ([5]) which, by defining p — pc = ua~, can be rewritten as 

u +-Pu + = 0. (29) 

To derive this relation we have used ^7 + 5 ( 7 )^^ = where P is the total pressure of the 

Universe. Using ~ —P 4 , ~ —V{pc), and the relation 

U + Ap ~ U (Pc) + A{pc)p + ^ {V4,4>{pc) + A^4,4,{Pc)p) {p - pcf , (30) 

which holds near the extremum of the potential pc, we finally obtain 

u + (^,4>4,iPc) + A^ 4 , 4 ,{pc)poa~^ - (pc)^ u = 0. (31) 


Setting A = 1 , this reduces to the result obtained in [5] and [TF]. Now in the regime 7 ^ ~ 0 , we 
approximate a(t) in ([3T]1 by (1271) . Changing the variables x = cosh(T) and u{x) = {x^ — l)iy{x), 
where r = ^ is a dimensionless time, (IHT]) reduces to the Legendre equation 

- A'Pp - + 1 ) - T^) = «■ ( 32 ) 

In terms of the scale factor, the parameter x is given by 


a; = ,/l _L ^<^0 a3 

r ^mo ' 


(33) 


where flmo = pm is the sum of the dark matter and baryonic matter densities. Also 

, AV'iPc) V{Pc) 1 

-v{v + l) = + 

3p0o P4>o 4 


= 7 


1 _ 4:A"{pc)pdo 
4 3A(^pc)pmO 


(34) 
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To get more physical intuition on ^ and v, we use 


V"{<Pc) + 


M4>c) 


'Pdc 0) 


(35) 


where pdc = Pd{o-c) is the value of dark matter density when the effective potential becomes convex. 
Using (l35]) . we rewrite ([M]) as 


-zy(i^ + l) = 


9Q 


<t>o 


V”{ct>c) \ _ 1 

Hi ) p<po 4’ 


= 


( 4ag \ 1 


\9nmoJ V Hi 


(36) 


// ^ 

with ^ ~ where is the mass of (f), e.g. for the potential (fTSll . = V2p which is the 


^0 


mass of excitations around the true minimum after the symmetry breaking. General real solution of 
the scalar field is given by [19] 


(37) 


</> = </>c + (^) - l)-3 {CiP-^{x) + C2e--^Q(;(x)) . 

V i ‘mO / 


Cl and C 2 are two real constants and Py and Qu are associated Legendre functions of the first and 
second kind respectively. Near w = —1, where (fP' <C V, we have 


102 _y(^) 

W 4 , = -TT— ^ -1 + 


i 02 + U( 0 ) 

From f we derive 

y/^mO 

Finally the EoS parameter of dark energy is obtained as 

A{p, V, x) 


V{^c)' 


(38) 




w^ = -1 + 


(x^ — 1)2 


(39) 


(40) 


with 


Aip,u,x)= [Cl + -)xP^ >^{x) - (u + p + l)P^^^{x) j + 
Cae--^ ((^2 + \)xQ^Ax) + {p - v - l)Q'',^\x)^ ]\ 


(41) 
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Dimensionless constants Ci are given by Cj = y 4 i/(^° )q° ^ Ci. The slow roll condition 02 <c U is valid 
when Wfj) = —1 + ~ _i. in terms of A (seeHDj), this can be rewritten as 


A{p,v,x) <C (x2 - 1 ) 2 . 


(42) 


This condition puts a constraint on the values of p and v in terms of the scale factor given by (|33l) 
in the slow roll regime. This constraint, as we will see, depends also on the value chosen for ac- To 
see wether a chosen p and u satisfies (|42F one can employ numerical computation. 

Since Ci are not independent, the boundary condition 


w, 


— Cic) — 1 , 


(43) 
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gives A{fi, u, ttc) = 0 , leading to 


^ (-/X - - 1) P^J^iiXc) + {V + 3/2) XcPu ^jxc) 

^ {fjL-u - l)Q'^_^_^{Xc) + + 3/2) XcQu{Xc) ^ 

where Xc = x{ac)- One may introduce another relation between the parameters of the model 

»,(a = l) = -l+(^)0(^,.,y^JT^)=«.„, (45) 


where wq is the EoS parameter of dark energy at the present time. Therefore if the parameters of 
the model are chosen according to (f36]l . dM]), ([H]) and (US]), the EoS parameter of the quintessence 
is determined by Wfp in (HHl) . 

To obtain more specihc relations for the parameters of the model, we may use astrophysical data. If 
based on Planck 2015 data [T], one chooses ~ 0.685, flmo — 0.315, and approximates V{4>c) ~ Pcpo 
(note that has the same order as pdo) , which is valid when <C P (or when w^p ~ —1), ([361) 
becomes 


P- 


= 7 + 1.41 
4 


i^(i^ + l) = 7-0.65 

4 


V'icPc 

-”o 

P"(0c 

Hi 


(46) 


Similarly ()^5]l reduces to 

- 1 + 0.312^(;u, z/, 1.78) = Wo, (47) 

where wq = —1.006 ± 0.045 with 95% C.L. [T]. 

Therefore, to study the evolution of wp, one needs to specify ^ and Oc. As an illustration let 
us take a model with ^ = — 10 “^, e.g. for (fTSIl this is equivalent to p = 0.0177o- We assume that 

the symmetry breaking occurred at Oc = 5 , and take wq = —0.99 [Tj. The EoS parameter of dark 
energy is plotted numerically in terms of x in figured] As it is clear from this figure < 0.01 so 


(+-l)a ¥(+) ^ 


X 



Figure 1: The EoS parameter of dark energy in thawing model in terms of x for ^ = —10 ^ and = h. 

Also, The EoS parameter of dark energy is plotted numerically in terms of the redshift z = ^ — 1 
in figure |21 
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Figure 2: The EoS parameter of dark energy in thawing model in terms of the redshift z for ^ —10 ^ and 

-”o 

flc = 2- 


5 Conclusions 

The thawing model has been considered in the literature previously [5] without considering any 
interaction between dark sectors. The evolution of the quintessence near the maximum of its own 
potential (not the effective potential) was also discussed in the absence of interactions in [8] as 
the hilltop quintessence. However it was not emphasized why the quintessence is near its potential 
maximum in the present epoch. 

In this paper, we considered the coupled quintessence in a thawing model and attributed the 
deviation of EoS parameter of dark energy from = — 1 and the dynamics of dark energy to the 
spontaneous Z 2 symmetry breaking of the effective potential. Indeed the beginning of the motion of the 
scalar field in the thawing model was related to the evolution of dark matter. We studied conditions 
required for the validity of this model for general potentials and couplings. We also investigated 
whether the scalar field can reside at the minimum of the effective potential and the stability of the 
model was discussed. 

Finally, in the last section, which includes our main results, we used the approximation W(p — — 1, 
and obtained analytical solutions for the quintessence and its EoS parameter near the extremum of 
the potential in terms of known functions. As far as possible we determined the parameters of the 
model in terms of astrophysical data. 

As an observational test and as an outlook, using equations ([57)1 . (1551) and (155)1 and the Friedmann 
equation one can try to obtain an expression for the Hubble parameter in terms of the redshift 
parameter 2 . By deriving H{z), one can obtain the luminosity distance redshift relation and compare 
it with type la supernova data [20] . 
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